High cumulants in the counting statistics measured for a quantum dot 
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Abstract 

We report on measurements of single electron tunneling through a quantum dot using a quantum point contact as 
non-invasive charge detector with fast time response. We elaborate on the unambiguous identification of individual 
tunneling events and determine the distribution of transferred charges, the so-called full counting statistics. We discuss 
our data analysis, including the error estimates of the measurement, and show that the quality of our experimental 
results is sufficiently high to extract cumulants of the distribution up to the 20th order for short times. 
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The progress in time-resolved charge sensing using 
quantum point contacts has made the detection of in- 
dividual electron tunneling events possible |Q1 0, 
Such experimental techniques enable measurements of 
the distribution of transferred charges in a quantum dot 
system, the so-called full counting statistics JUS 1(1. 
Measurements of the full counting statistics have been 
used to study the first few higher order cumulants of 
the distribution [7]. The first cumulant is the mean of 
the number of transferred electrons n, the second is the 
variance, and the third is the skewness. Recently, the 
general behavior of cumulants of very high orders was 
investigated experimentally as well as theoretically Jit] . 

Here, we report on real time single electron counting 
with a large bandwidth detector. This enables us to mea- 
sure the full counting statistics for single electron trans- 
port through a quantum dot with high precision. From 
the distribution of tunneling events the cumulants de- 
scribing the statistical properties of the system are ex- 
tracted. The cumulants of the transport statistics are 
widely studied in theory, and here we discuss the ex- 
perimental precision at which they can be measured. 

Our device is based on a GaAs/AlGaAs heterostruc- 
ture containing a two-dimensional electron system 
(2DES) 34 nm below the surface. The electron den- 
sity is p — 4.59 • 10 15 irT 2 , and the mobility is yiz = 
64.3 m 2 /Vs. We have used an atomic force microscope 
(AFM) to define the quantum dot (QD) and the quantum 
point contact (QPC) structure by local anodic oxidation 
(LAO) on the surface fiHU; the 2DES below the ox- 
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idized surface is depleted and insulating areas can be 
written. 

A schematic view of our device is presented in 
Fig. [T](a). The bright lines depict the insulating oxide 
barriers written by the AFM. The QPC (right area) is 
separated from the QD structure (left area) by an insu- 
lating line. The QPC can be electrically tuned using an 
in-plane gate by applying the potential Vc3- The cur- 
rent through the QPC is amplified with a 100 kHz cur- 
rent amplifier and detected in a time-resolved manner. 
A small bias is chosen to avoid back-action on the QD 
1 1 ill . The QD is coupled to source and drain electrodes 
via two tunneling barriers, which can be separately con- 
trolled with gate voltages Vgi and Vci- These gates are 
also used to tune the number of electrons on the QD. 
A sufficiently large bias V$d ensures that transport is 
unidirectional. The sample is placed in a dilution refrig- 
erator with temperatures in the mK-regime. 

The quantum dot is operated in the Coulomb block- 
ade regime close to a charge degeneracy point, where 
a single electron at the time can enter the QD from the 
source electrode and leave via the drain. Consequently, 
the electron number on the QD changes back and forth 
between and 1 additional electrons as electrons tun- 
nel through the system. This leads to a corresponding 
change of the electrostatic potential at the position of the 
QPC. The QPC is tuned to a working point on the edge 
of the first conductance step, where the current through 
the QPC is very sensitive to the presence of localized 
electrons on the QD. With 1 additional electron on the 
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Figure 1 : Operating principle for measuring single electrons passing 
through a QD using a QPC as detector, (a) Schematic picture of the 
experimental setup. Current through QPC and QD are measured in 
a two circuit setup. The QPC current is measured with a 100 kHz 
bandwidth, (b) Dots: Time segment of the raw measured QPC current 
signal. Line: Resulting signal after analysis of the raw data. Tunneling 
times Ts and td deduced from this analysis are marked by dotted lines. 
Histogram: Distribution of currents from the time trace. From the two 
peaks in the distribution, two bands are extracted corresponding to 
(upper band) and 1 (lower band) additional electrons on the QD. 



QD, the current through the QPC is suppressed, and the 
suppression is only lifted as the electron leaves the QD 
via the drain electrode. 

The time-resolved QPC current consequently shows 
a random telegraph signal, as shown in Fig.Q](b). The 
individual data points (dots) correspond to an integra- 
tion time of 10 /vs. Most of the data is distributed in 
two bands, or current levels. The histogram of the time 
trace on the right shows the two levels of the QPC cur- 
rent. The lower (upper) state corresponds to the QD 
with 1 (0) additional electrons. The time intervals that 
the QPC current stays in either of the two bands directly 
yield the waiting times for an electron to tunnel on (js ) 
or off (td) the QD. The tunneling events themselves are 
marked by the transitions between the two bands. 

In order to access the transport statistics it is neces- 
sary to extract the number of tunneling events from the 
raw data. This is done using an averaging algorithm. 
The key principle is to identify the two bands corre- 
sponding to the different charge states of the QD. To 
this end we use the histogram of the QPC current. The 
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Figure 2: (a) Distribution compiled from 853,181 tunneling events. 
The experimental result (bars) are compared with a Poisson distribu- 
tion (dashed line), a Gauss distribution (dotted line), and theoretical 
calculations corresponding to the experimental setup (solid line), (b) 
Relative deviation 6N between the experimentally determined distri- 
bution and the theoretical distributions in panel (a). 



two peaks are identified and the boundaries of the two 
bands are located. Once the boundaries of the upper and 
lower band are known, tunneling events can be identi- 
fied as transitions between the two current levels. We 
identify a transition by checking each measured current 
point for being part of one of the two bands. As long as 
the system stays in one state, the current level is aver- 
aged. In order to avoid a shifting average due to external 
noise sources, we neglect points that are not part of one 
of the two bands. If a point is located in another band 
than the previous point, the charge on the quantum dot 
has changed, and a tunneling event is registered. 

The band selective averaging approach features a 
number of advantages. Most importantly, it is more 
robust to external disturbances of the measured current 
than the alternative method of edge detection directly at 
the flank of a transition. As only points inside the two 
bands are taken into account, noise spikes are not affect- 
ing the detection of a transition. The current averaging 
itself gives an additional verification, as a shift of the av- 
erage between two tunneling events can be detected. In 
this way, unintended changes of the detector are directly 
observable. Also data points outside the bands can be 
used to identify the power of external noise added to the 
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Figure 3: High order cumulants (a) Normalized cumulants ((n"'})/((n}) for m = 2-5. Curaulants are shown for times fo up to 3 ms. (b) Cumulants 
for m = 10, 15, 20 for times fo up to 750 /.is. (c) Error estimates calculated as the difference between cumulants obtained from the first and second 
half of the full time trace. The smoothened absolute error is shown for direct comparison with the measured cumulants. The error estimates for 
the measured cumulants (a, b) are shown with thick lines, while the thin lines correspond to the error estimates for other cumulants up to order 20. 
Where the solid lines cross the dashed line the corresponding estimated error is ~ 10%. 



detector current. This allows us to identify time seg- 
ments with too high noise level and it also improves the 
effective bandwidth of the detector. 

We now turn to the statistics of the number of charges 
transferred through the QD. The exact times of tunnel- 
ing events are stochastic, and the statistical properties 
of the transfer process are best captured by the distribu- 
tion of the number n of electron counts within a fixed 
time segment length to- An example of such a dis- 
tribution, compiled from more than 800,000 tunneling 
events, is shown in Fig. [2] (a). Here, fo = 20 ms, with 
a corresponding mean value ((«)) of 17.13 electrons per 
time segment. The measured distribution (bars) is com- 
pared to a Poisson distribution (dashed line) with the 
same mean as the experimental data, a Gauss distri- 
bution (dotted line) with the same mean and variance 
as the experimental data, and theoretical calculations 
of the distribution based on a model of the experimen- 
tal setup (solid line) [8]. The relative deviation 6N of 
these curves compared to the measured distribution is 
depicted in Fig.[2](b). 

The Poisson distribution deviates significantly from 
the measured distribution due to its larger width. 
The sub-poissonian width of the measured distribution 
demonstrates noise suppression due to Coulomb block- 
ade on the QD [ 12]. The Fano factor, i.e., the variance 



normalized with respect to the mean (which is 1 for a 
Poisson process), is reduced to nearly 0.5 due to the 
almost symmetric tunneling rates. For the Gauss dis- 
tribution, the variance enters as an additional parame- 
ter, and it clearly fits better than the Poisson distribu- 
tion. Still, however, noticeable deviations are seen in 
the tails of the distribution. The best agreement with ex- 
periment, in particular in the tails of the distribution, is 
given by the theoretical curve derived from a two-state 
rate model, taking into account the finite detector band- 
width of about 40 kHz. The rates used for this calcula- 
tion were determined from the waiting time distribution 
as described in Ref. J^]. 

Next, we examine the cumulants of the distribution. 
As the distribution depends on the time segment length 
to, also the cumulants are time dependent. The cumu- 
lants of order 2 to 5, normalized with respect to the 
mean, are shown in Fig. [3] (a). All normalized cu- 
mulants start from a value of 1 at very short times, 
where the transfer process is nearly poissonian. For 
a Gauss distribution, only the first and second cumu- 
lants are non-zero. We observe a clear non-zero value 
for the third cumulant, showing that the distribution is 
not Gaussian, in agreement with the observed devia- 
tions in Fig. [2] For longer times the normalized cu- 
mulants reach a constant value, referred to as the long- 



time limit. While the second and third cumulants show a 
monotonic transition between 1 and the long-time limit, 
the fourth and fifth cumulants have a minimum at short 
times. Additionally, the fifth cumulant has a slight max- 
imum, following the minimum, showing some indica- 
tions of an oscillatory behavior. This phenomenon is 
further enhanced with increasing orders of the cumu- 
lants as seen in Fig.[3](b), showing cumulants of order 
m = 10, in — 15 and m = 20 that clearly oscillate with 
time. Not only does the number of oscillations increase, 
but the amplitude also grows dramatically. For m = 10, 
the amplitude reaches a value of 16, for in — 15 the am- 
plitude is around 10,000, and for in = 20 the amplitude 
exceeds 40,000,000. As a theoretical analysis has re- 
cently revealed, the magnitude of the oscillations grows 
factorially with the cumulant order 01 . More surpris- 
ingly perhaps, the oscillatory behavior of the cumulants 
is to be expected in a large class of nontrivial distribu- 
tions as functions of almost any parameter 0. 

The errors on the measurement increases with the or- 
der of the cumulants and the time segment length to. 
This is supported by the simple error estimates in Fig. 
[3(c), where an absolute value for the estimated error is 
shown as functions of the cumulant order and the time 
segment length. The error estimates are obtained by 
splitting the full data set into two pieces and calculating 
the cumulants for both halves individually. The absolute 
value of the difference is then used as a simple estimate 
of the statistical error. The result is smoothened using 
a box averaging. For increasing time an approximately 
exponential increase of the error is observed for each 
cumulant. A dotted line indicates the time at which the 
estimated error reaches about 10% of the absolute value 
of a given cumulant. In this experiment we were able 
to measure the long-time limit of cumulants up to order 
m = 7, which was reached after 2-3 ms. 
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In conclusion, we have performed measurements of 
the full counting statistics of electrons passing through 
a quantum dot. We have described our detection al- 
gorithm, shown experimental and theoretical results 
for the distribution functions, and presented measure- 
ments of the time evolution of the high-order cumu- 
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